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Fluctuations and all-in—all-out ordering in dipole-octupole Nd 2 Zr 207 
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We present an experimental study of the pyrochlore coumpound Nd 2 Zr 207 by means of neutron 
scattering and magnetization measurements down to 90 mK. The Nd^^ magnetic moments exhibit 
a strong local (111) Ising anisotropy together with a dipolar-octupolar nature, different from the 
standard Kramers-doublet studied so far. We show that, despite the positive Curie-Weiss tempera¬ 
ture, Nd 2 Zr 207 undergoes a transition around 285 mK towards an all-in—all-out antiferromagnetic 
state. We establish the {FI, T) phase diagram in the three directions of the applied field and reveal 
a metamagnetic transition around 0.1 T. The strongly reduced ordered magnetic moment as well 
as the unexpected shape of the magnetization curves demonstrate that Nd 2 Zr 207 is not a stan¬ 
dard Ising antiferromagnet. We propose that the peculiar nature of the Nd doublet combined with 
competing interactions explain these findings. 


Understanding, characterizing and classifying novel 
states of matter is one of the main issues of research 
in frustrated systems mm- Much of the recent activ¬ 
ity in this field has focused on 4/ pyrochlores. Here, the 
frustration arises from the lattice, built from tetrahedra 
connected by their corners. The situation where the mag¬ 
netic moments on these vertices are coupled via a ferro¬ 
magnetic nearest neighbors interaction and constrained 
by the crystal electric field (GEE) anisotropy to lie along 
local (111) axes, has retained a lot of attention, as it leads 
classically to a macroscopically degenerated ground state, 
where each tetrahedron has two spins pointing in and two 
spins pointing out, called “spin ice” mm- The quest for 
the quantum variant of spin ice (QSl) has attracted fur¬ 
ther strong interest. A wealth of fascinating phenomena, 
like the emergence of gauge fields along with unconven¬ 
tional excitations are expected [S] . Those systems yet re¬ 
main to be explored from an experimental point of view. 
In this search, systems with positive Curie-Weiss tem¬ 
perature dew indicative of resulting ferromagnetic inter¬ 
actions along with strong effective Ising-like anisotropy 
are a priori possible candidates. In this letter, we fo¬ 
cus on Nd2Zr207 which falls into this class, as shown by 
magnetization EHU and CEE studies mm- 

In addition, Nd2Zr207 stands out through an interest¬ 
ing property. In many of the 4/ pyrochlores, the crystal 
field and spin orbit coupling lead to Kramers ground dou¬ 
blets, well separated from the excited states, so that the 
low energy physics can be described by a pseudo spin 
1/2 [To]. In most of them, this spin 1/2 transforms as a 
magnetic dipole under space group operations, just like 
a true spin 1/2 moment. This is for instance the case 
in Er2Ti207 HH [T2], Er2Sn207 [T^ or Yb9Ti907 [TU] . 
However, there exists Kramers doublets called “dipolar- 
octupolar” (DO) doublets, where, in a pseudo spin 1/2 
representation, two components transform like a mag¬ 
netic dipole, while the other transforms like the com¬ 


ponent of the magnetic octupole tensor [T3]. This is 
the case in Nd2Zr207 |S] (but also in Nd2lr207 [15] 
and Dy2Ti207 HU). The theoretical phase diagram of 
such systems is highly interesting m- besides ordered 
phases like the all-in—all-out structure (AIAO) where 
dipoles point along the local axes toward or away from 
the tetrahedron centers, an antiferro-octupolar phase ap¬ 
pears along with QSI phases. 

In this letter, we first confirm the strong Ising charac¬ 
ter of the Nd^+ ion in Nd2Zr207 and the DO nature of 
its ground state doublet. Combining neutron and magne¬ 
tization measurements, we show that despite a positive 
(^cw, Nd2Zr207 orders in an antiferromagnetic AIAO 
structure below Tn = 285 mK with a strongly reduced 
magnetic moment. We then establish the {H,T) phase 
diagram: a metamagnetic-like transition occurs around 
0.1 T, and a conventional field induced ordered state is re¬ 
covered at about 0.2 T. This can be reproduced by mean 
field calculations with an effective antiferromagnetic cou¬ 
pling that however contradicts 0cw > 0. The reduced 
ordered magnetic moment and the peculiar shape of the 
magnetization versus field curves can be accounted for by 
considering local fluctuations induced by octupolar cor¬ 
relations. We propose that, despite the presence of fer¬ 
romagnetic interactions inferred from positive 0cw that 
should stabilize the spin ice state, the peculiar nature of 
the ground state doublet combined with competing in¬ 
teractions induce the AIAO structure. 

Magnetization and ac susceptibility measurements 
were performed in two SQUID magnetometers equipped 
with dilution refrigerators developed at the Institut Neel 
[18] . The experiments were carried out on a powder 
sample and a parallelepiped single crystal of Nd2Zr207 
grown by the floating zone technique (2.35 x 1.72 x 1.36 
mm^) [19] down to 90 mK. The magnetic field was ap¬ 
plied along the three cubic high symmetry directions 
[100], [110] and [111]. Powder neutron diffraction ex- 


2 


B 20 

Bio 

B 43 

Beo 

Bqs 

Bee 

9\\ 

91- 

-190 

5910 

110 

2810 

10 

-890 

4.3 

0 


A 

-9/2 - 

7/2 -5/2 -3/2 ±1/2 

3/2 

5/2 7/2 

9/2 

It) - 

-0.878 

0 0 

-0.05 

0 

0.476 

0 0 

0.009 

U) - 

■0.009 

0 

0 

0.476 

0 

0.05 

0 0 

-0.878 


TABLE I. Stevens coefficients (in K), g-factors and ground 
state wave-functions of Nd®^ in Nd 2 Zr 207 from inelastic neu¬ 
tron scattering analysis. 


periments have been conducted on DIB (CNRS CRG - 
ILL) using a wavelength A = 2.52 A. We repeated the 
same experiment on a better quality sample (the one 
used for macroscopic measurements) on 4F2 (LLB), using 
A = 2.36 A. Single crystal diffraction experiments were 
carried out on D23 (CEA CRG - ILL) using A = 1.28 
and 2.4 A. Finally, CEF excitations were measured on 
the high quality powder sample by means of inelastic 
neutron scattering conducted on the thermal triple axis 
spectrometer 2T (LLB). 

First, we focus on the DO nature of the ground dou¬ 
blet in Nd 2 Zr 207 and study the CEF excitations. As 
shown in Figure [^a), the spin-spin correlation function 
S{Q,uj) exhibits two clear peaks at 21 and 33 meV. The 
intensity of both excitations tends to slightly decrease 
as the temperature increases and remains roughtly con¬ 
stant at larger wavevector Q. These two characteristics 
are evidence for the magnetic origin of both modes, that 
are thus attributed to two CEF transitions. A better 
energy resolution reveals the presence of an additional 
single mode at 10 meV. Its temperature and Q evolu¬ 
tion, strong increase with increasing temperature or in¬ 
creasing Q, are typical of a phonon mode. The energy 
of these two CEF modes match the theoretical analy¬ 
sis based on macroscopic measurements, that takes into 
account the mixing of the “^ 1^/2 and ‘^Iii /2 multiplets of 
the 4/^ Nd^+ ions electronic configuration [5]. Note that 
the theory predicts two distinct modes around 32 meV 
which could not be resolved in the experiment [20] • To 
explicitely show the DO character of the doublet, and 
determine wavefunctions easier to handle, we consider 
the CEF Hamiltonian within the ground multiplet ‘^Iq /2 
only: = Y.m,n B-nmOnm where the Onm are the 

Stevens operators (J = 9/2, gj = 8/11) [21]. The Bnm 
coefficients are refined by fitting the neutron data (en¬ 
ergy levels and intensities) along with the value of the 
local magnetic moment at low temperature 2.3 /is (see 
Figure [^a) and Table 0- This procedure gives modes at 
22, 33.5, 34.5 and 109 meV yielding a picture of very 
strong Ising Nd^+ moments with g_L = 0 and g// = 4.3 
in agreement with Ref. [8] . The coefficients of the ground 
doublet wave-functions \cj =t,4.) (see Table[l| are similar 
to the ones found in Nd 2 lr 207 m and typical of a DO 
doublet [11] . 
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FIG. 1. (a) Inelastic neutron scattering spectra taken at 

several temperatures. The solid line shows the fit with the 
coefficients of Table [^ (b) Powder diffractogram at 4 K. The 
line is a fit to the FdZm structure. Inset: Scheme of the AIAO 
state, (c) Zoom on the low angle part at 150 mK and 4 K 
showing the magnetic contribution. The lines are combined 
refinements of the crystallographic and magnetic structure. 


Our magnetization and ac susceptibility measurements 
reveal a transition at low temperature, as expected from 
specific heat measurements [S]. The susceptibility ex¬ 
hibits a peak both in the real and imaginary parts at a 
temperature of 410 mK for the powder sample, consistent 
with ref. [6] (not shown), and at a lower temperature of 
285 mK for the single crystal (see Figure]^. Despite 
0CW > 0 (see inset of Figure [^, the shape of the y' 
peak, as well as its amplitude are characteristic of an 
antiferromagnetic transition. The position of the peak 
does not depend on frequency (in the measured range 
/ = 0.11 — 570 Hz). However, when the frequency in¬ 
creases, the peak of y" lies on top of a background below 
4 K which increases with decreasing temperature, simi¬ 
larly to what was observed in Tb 2 Ti 207 [22] (not shown). 
This observation might be the signature of slow fluctua¬ 
tions, not affected by the antiferromagnetic transition. 

This antiferromagnetic transition is further confirmed 
by neutron diffraction measurements. At 4 K (see Fig¬ 
ure [^b)), the FULLPROF refinement [23] in the Fd3m 
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FIG. 2. ac susceptibility Xi and vs temperature T along 
[110] with iFac = 2.7 Oe and / = 5.7 Hz, corrected for demag¬ 
netization effects. Inset: l/xiVsT. The line is a fit between 1 
and 4 K to the Curie-Weiss law l/y = 2.05x 10"^ —1.05x 10® T, 
giving Ocw = 195 (±15) mK and an effective moment = 
2.45 (±0.02) HB- 


space group yields a lattice parameter a = 10.66 A and 
the 02 oxygen position x = 0.335 (with an RF factor of 
2). A magnetic contribution rises below Tjq and appears 
on the nuclear Bragg peaks, indicating that the propa¬ 
gation vector of the magnetic structure is k = ( 0 , 0 , 0 ) 
(see Figure [^c)). The symmetry analysis shows that a 
number of magnetic configurations, corresponding to the 
different irreducible representations ± 3^5 7 _g are in princi¬ 
ple possible candidates [53]. Fa corresponds to the AIAO 
structure (see inset of Figure [^b)) predicted for multi¬ 
axis Ising pyrochlore antiferromagnets [25], ±5 to the 
structure of the antiferromagnetic XY pyrochlore magnet 
ErgTigOy [55], ±7 to the structure of dipolar pyrochlore 
magnets such as Gd 2 Sn 207 [53 [5H], and Fg to the splayed 
ferromagnetic Yb 2 Sn 207 structure [55] [55]. Our neu¬ 
tron data are clearly compatible with Fg since it is the 
sole configuration that leads to a vanishing contribution 
of the (111) and (002) (plus symmetry related) Bragg 
peaks, as observed in the experiment. The fullprof re¬ 
finement of the diffractogram yields a weak magnetic mo¬ 
ment m = 0.8±0.05 fiB indicating large fluctuations even 
at temperatures of about T^/A. Diffraction measure¬ 
ments performed on a single crystal confirm this analysis 
and lead to the same conclusion. Note that this AIAO or¬ 
dering, although classically expected in Ising pyrochlore 
antiferromagnets, has rarely been reported. It has been 
observed in compounds such as FeFa [30], Cd20s2Sb7 
m or Na 3 Co(C 03 ) 2 Cl [35], but in these cases, it does 
not result from Ising anisotropy. The ordering of the 
rare-earth moments in pyrochlore iridates |33[ 134] may 
be placed in this category, but this ordering is induced 
by the AIAO state of the iridium magnetic moments and, 
a priori, not to the magnetic interactions between rare- 
earth moments. In that sense, the observation of the 
AIAO state itself makes Nd 2 Zr 207 a peculiar case in the 
insulating pyrochlore compounds. 



FIG. 3. Magnetization M vs internal field Hi measured along 
the three main directions at 90 mK. Inset: Detailed view of 
the low field data. The field was swept up from negative to 
positive values. 


Aiming at a full characterization of this system, we 
proceed with the determination of the {H, T) phase di¬ 
agram by performing M ys H measurements down to 
90 mK. As previously reported [8], the saturated magne¬ 
tization depends on the field direction (see Figure]^. The 
obtained values are, within 10%, the ones expected for 
magnetic moments with a multi-axis Ising anisotropy |35j 
and correspond to a total magnetic moment m « 2.3 /is, 
consistent with the effective moment of the Curie-Weiss 
law (see inset of Figure]^. This shows that, although the 
ordered moment is reduced in the AIAO phase, the total 
magnetic moment is recovered in high magnetic fields, 
and remains with a strong Ising character. At low field, 
the slope of the magnetization curves is not zero (see inset 
of Figure [^, contrary to what is expected for a strongly 
anisotropic antiferromagnet. This confirms that a part of 
the magnetic moment is still fluctuating at low tempera¬ 
ture. An inflexion point is present in these curves around 
0.1 T, which can be attributed to a metamagnetic process 
towards the field induced ordered phase. It is observed 
in the three directions of the applied field, but is more 
pronounced along the [111] direction. The temperature 
dependence of this metamagnetic process can be followed 
by plotting the derivative dM/dHi of the magnetization 
curve (see inset of Figure]^. Assuming that the maxi¬ 
mum of dM/dHi corresponds to the metamagnetic field 
value, we plot the phase diagram as a function of tem¬ 
perature (see Figure]^. As expected the dM/dHi peak 
disappears at the transition temperature. The obtained 
phase diagram is almost isotropic, and the data can be 
fitted to a T — Tjv law which gives an exponent slightly 
larger than 1/2. 

We now compare these experimental results with 
mean-field calculations. To this end, we consider a model 
taking into account 'HQgp and a unique coupling act- 
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FIG. 4. Phase diagram {H, T) obtained from the maximum of 
dM/dHi vs Hi. Inset: dM/dHi vs Hi in the [111] direction. 


ing between nearest neighbor local z components of the 
Nd moments. An antiferromagnetic ~ —0.005 K al¬ 
lows to reproduce the Neel temperature of about 300 
mK and the AIAO structure. This leads to a {H, T) 
phase diagram similar to the experimental one, with a 
metamagnetic-like transition from the AIAO state to the 
field induced ordered state (3-in—1-out or 2-in—2-out 
state depending on the direction of the field). Neverthe¬ 
less, this mean-field model fails to reproduce the positive 
0 CW as well as the small amplitude of the ordered mo¬ 
ment. In addition, the calculated field induced transition 
is abrupt, corresponding to spin-flips in the whole mag¬ 
netic structure, in contrast with the smooth measured 
curves (see dashed lines in Figurea)). 

These results support the existence of strong fluctu¬ 
ations which prevent a total ordering of the moment. 
Further, from the coefficients of the DO wavefunction, 
(t I J| i) = 0, so the dipolar exchange alone cannot in¬ 
duce any on-site fluctuation from one element of the dou¬ 
blet to the other. This suggests that another coupling is 
at play. Considering the DO nature of the ground wave- 
function, such fluctuations arise by virtue of the octupo- 
lar operator T = i{J'^— J~J~J~) [T7j. Indeed, 
using Table |Tj one has (f |T| i) 7 ^ 0. As a result, an oc- 
tupolar field induces a mixing term in the Hamiltonian 
V = X]o'=t 4 which tends to bind the two el¬ 

ements of the doublet to form new states with reduced 
moments. In this picture, applying a magnetic field com¬ 
petes with V and progressively unbinds | f) and | |). The 
effect of such term is consistent with the experimental ob¬ 
servations, since, as shown in Fig. it does reduce the 
ordered moment, and makes the magnetization curves 
and field induced transition smoother. However, the cal¬ 
culation of the susceptibility with this value of and 
taking into account V still cannot reproduce the sign of 
(^cw- 

This discrepancy might be understood in the con¬ 
text of the pseudo spin 1/2 theory developed recently 
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FIG. 5. (a) Calculated M vs H at T = 100 mK with ~ 
—0.005 K for the field applied along the [100] (green) and 
[111] (red) directions. Full (resp. dashed) lines are calculated 
with U = 0.2 K (resp. U = 0). (b) Calculated Phase diagram 
with H II [100] superimposed with M. ~ —0.005 K and 
U = 0.2 K. The dashed line shows the transition. 


for DO doublets m- It is shown that despite a posi¬ 
tive (corresponding to a ferromagnetic-like interaction 
in that local frame convention), the antiferromagnetic 
AIAO ground state is stabilized by a negative transverse 
term in a large region of the theoretical phase dia¬ 
gram. But due to the peculiar nature of the ground state 
doublet, only contributes to the Curie-Weiss tempera¬ 
ture, leading to a positive 0cw- Further analysis of spin 
waves would help to determine these coupling constants 
and serve as a test of this hypothesis. 

It is worth noting that this contradiction might also be 
removed in the context of the recently proposed fragmen¬ 
tation theory |36j . In this model, for a certain range of 
microscopic parameters, the classical spin-ice phase can 
fragment into an AIAO antiferromagnetic phase with a 
reduced moment, and a Coulomb phase characteristic of 
the spin-ice correlations. Such a scenario could be con¬ 
firmed by diffuse scattering measurements, showing the 
coexistence of Bragg peaks and pinch points. 

In summary, our inelastic neutron scattering results 
confirm the Ising and DO nature of the ground state dou¬ 
blet of the Nd^+ ion. Despite a Curie-Weiss temperature 
indicative of ferromagnetic interactions, the existence of 
an AIAO ordering is established. Strong fluctuations, 
unexpected in an Ising magnet, reduce however the mag¬ 
netic moment. We also provide evidence for a metam- 
agnetic transition towards a field-induced ordered state 
below the Neel temperature. More generally, this study 
points out how quantum fluctuations are reintroduced in 
Ising systems by multipolar correlations. In Nd 2 Zr 207 , 
the DO nature of the doublet naturally calls for octupolar 
correlations. Further theoretical and experimental work 
is needed especially to determine the nature of magnetic 
excitations stemming from this peculiar ground state. 
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